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Abstract. We study the traveling wave solutions to a reaction diffusion sys- 
tem modeling the public goods game with altruistic behaviors. The existence 
of the waves is derived through monotone iteration of a pair of classical upper- 
and lower solutions. The waves are shown to be unique and strictly monotonic. 
A similar KPP wave like asymptotic behaviors are obtained by comparison 
principle and exponential dichotomy. The stability of the traveling waves with 
non-critical speed is investigated by spectral analysis in the weighted Banach 
spaces. 



1. Introduction 

Due to their widespread applications in biology, physics and chemistry, there 
has been an increasing interest in the traveling wave solutions to reaction diffusion 
systems. A fruitful methods have been developed in deriving the traveling wave 
solutions, among which the monotone iteration method is proved to be rather ef- 
fective. Such method reduces the existence problem to that of an ordered pair of 
upper and lower solutions. However, it is not an easy task to construct the upper 
and lower solution pairs in the classical sense, therefore, [4, 19, 25, 26J turn to search 
for the generalized upper and lower solution pairs and have achieved the desired 
results. Other methods such as geometric singular perturbation [213, phase plane 
analysis [H [U US [THl HH HS1 HH HH Ell [2H H3l ISl] and the fixed point Theorem 
[19] are also successfully applied to various reaction diffusion systems to obtain the 
wave solutions. 

One interesting yet elusive question is the long term behavior of the traveling 
wave solutions. In fact, even for the relatively well studied Lotka Volterra systems 
less results (see pQ, [T2]) are known on the asymptotic decay (growth) rates of 
the waves and their stability. This question seems to be related to environmental 
selection in a given ecosystem. Roughly speaking, traveling wave is formed by 
intersecting the unstable subspace from one equilibrium and stable subspace from 
another. The stable and unstable spaces are generally multi-dimensional, therefore, 
the traveling wave selects an unstable direction from the unstable equilibrium to 
emanate, and then a stable direction from the stable equilibrium to enter. The 
selection affects the asymptotic behaviors and stability of the wave solutions. 
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We try to understand those problems through the study of the traveling wave 
solutions of the following reaction diffusion system: 



(1.1) 
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dt dx' 2 
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where u = u(x,t), v — u(x,t) are two competing populations and n, r2, a are 
positive constants. The system is a continuous spatial-temporal version of a public 
goods games (24 ]) and describes the interaction between the populations u(x,t) 
and v(x,t). In the game, the population u employs an altruistic strategy in order 
to ensure the survival of both populations; while the population v's strategy is to 
maximize its own gain. The function k{u) = ko + ku with k, kg > represents the 
public goods, contributed by the population u and shared by the population v. In 
this sense, we call population u the cooperators and population v the defectors. The 
cooperators have to pay certain penalty measured by au for such altruistic strategy. 
As the outcome of the game, the population v eventually wins while the population 
u loses. So two equilibrium states are of interest to us - the one dominated by u 
and the one dominated by v. The game is played by finding a passage between the 
two equilibrium states. 

After rescaling the density functions (see [9]), we may further assume that n = 
r'2 = fco = 1- The rescaled system has two unstable constant equilibrium states 
(0,0), (K*,0) with K* = l l^+ ak > and one asymptotically stable state (0,1). 
Throughout the rest of the paper we make the following assumption: 

H: < a, k < 1. 

The transformations u = K* — u, v = v change (|1.1[) into a local monotone type. 
Disregarding the initial conditions, the traveling wave solution to (jl . 1|) has the form 
(«(£), v(£)) T := (u(x + ct), v(x + ct)) T , c > and solves the following system 

u" - cv! - (K* -u){l ^-u+v 



(1.2) 



cv' + v(l — 



K ' 



l+k(K 



l+k{K* 



^t) =0, 



a) = 0, 



(u,v) T (~oo) = (0,0) T , (u,v) T (+^) = {K\l) T . 

where T represents the vector transposition and the prime the derivative with 
respect to £. 

In our derivation of the traveling wave solutions, the classical results on KPP 
(Kolmogorov-Petrovskii-Piskunov, |18j ) equation and the monotone structure of 
(|1.2[) play a central role. Noting that if u is replaced by zero then the second 
equation in (|1.2p is a KPP equation. The construction of the upper- and lower 
solution pairs is based on this simple observation. Cares also be taken on the 
asymptotic decay rates of the upper and lower solution pairs to ensure that they 
are ordered and have the same decay rate at negative infinity. The existence of the 
traveling wave solutions for system (|1.2| then follows from the monotone iteration 
scheme by [4] and [26j , and the asymptotics of the traveling wave will be obtained 
from the exponential dichotomy and comparison. 

The stability of the traveling wave solutions with non-critical speed is investi- 
gated by spectral analysis. Since the wave solutions are essentially unstable in the 
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unweighted spaces such as L°° or C un jf , similar to the stability issues of KPP waves 
[TU1 |2"Tl |2"TI 121)] , we first localize the initial perturbations to certain exponentially 
weighted Banach spaces. Then following the standard treatment of such problem 
as in [2] [17l [20] [21] , we show that the linearized operator around the wave solution, 
with noncritical wave speeds, does not have eigenvalues with positive or zero real 
part. Several new ideas are incorporated to locate the eigenvalues of the linearized 
operator. This, along with the fact that the linearized operator is sectorial leads to 
the nonlinear stability of the wave solutions. The stability of the wave solution im- 
plies that a sufficiently small initial disturbance will not change the outcome of the 
game. We remark that though the paper is written for a particular system the 
method presents here can be easily adapted to other reaction diffusion systems with 
monstable structure. We further remark that the stability of the waves with critical 
wave speed has a complete different nature and requires a different treatment, so 
we leave it open for now. 

The paper is organized as follows: in Section 2 we show the existence, uniqueness, 
monotonicity, and derive the asymptotic decay (growth) rates of the traveling wave 
solutions. In Section 3 we study the stability of the traveling wave solutions in 
weighted Banach spaces. 

2. Existence and Asymptotic Decay Rates 

In this section, we establish the existence of traveling wave solutions and derive 
their asymptotic decay rates through monotone iteration of a pair of smooth upper 
and lower solutions. The definition of the upper and lower solutions is standard. 
For the rest of the paper the inequality between two vectors is component- wise. 

Definition 1. A smooth function v(£ )) T , £ G R is an upper solution of (|1.2|) 

if its derivatives (u',v') T and (u",v") T are continuous on R, and (u, v) T satisfies 

f u" - cu' - (K* - u)(l - ^(7+;) - a) < 0, 
(2-1) \ 

1 «"-^' + «(i-tS^)<o, 

with the boundary conditions 

m (:)<->-( S) - (T)- 

A lower solution of (|1.2[) is defined in a similar way by reversing the inequalities 
in d5T]) and (|2~2|) . 

The construction of the smooth upper and lower solution pairs is based on the 
solutions of the KPP equation: 

( w" - cw' + f(w) = 0, 

(2.3) I 

[ w(— oo) = 0, w(+oo) = b. 

where / € C 2 ([0, b}) and / > on the open interval (0,6) with /(0) = f(b) = 0, 
f (0) = ai > and f(b) = -b x < 0. We first recall the following result ([21]): 

Lemma 2. Corresponding to every c > system (|2.3p has a unique (up to a 

translation of the origin) monotonically increasing traveling wave solution it>(£) for 
£ G R. The traveling wave solution w has the following asymptotic behaviors: 
For the wave solution with non-critical speed c > 1^fa\, we have 

(2.4) w (£,)=a w e C °* ^ * + o{e !L ^ lj: ^ ±i ) as £ -> -oo, 
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(2.5) w(g)=b-b w e 5 4 + o(e 5 6 ) as £ -> +oo, 

where a w and b w are positive constants. 

For the wave with critical speed c = we have 

HO = ( a c + dcO eV514 + °(£e^«) as £ -» -oo, 



where the constant d c is negative, b c is positive and a c eR. 

To construct the upper-solution for the system (| 1 - 2[) , we begin with the following 
form of KPP system: 

— v(l-v) = 0, 



(1 - k+ ak)(l + kK*(l - v)) 



(2.6) 

v(— oo) = 0, 5(+oo) = 1, 

where relating to f(v) = {1 ^ k+ak)[ i +kK « (1 -. d)) v{l - v) > for v G (0,1). 

/(0) = /(l) = 0, /'(0) = a > and /'(l) = -a/(l - fc + afc) < 0. According to 
LemmaJH for each fixed c > 2y/a, system (|2.6p has a unique (up to a translation of 
the origin) traveling wave solution v(£) satisfying the given boundary conditions. 
Define 

/ u(0 \ ( K*v(0 \ 
(2.7) = , (£K, 

V m J V m J 

then we have the following result, 

Lemma 3. For each fixed c > 2y / a7, (|2.7|) is a smooth upper solution for system 
(O). 



Proo/. On the boundary, one has (u,"u) T (-oo) = (0, 0) T , (u, v) T (+oo) = (K* , 1) T . 
As for the u component, we have 

, , w K*-K*v + v, 
u" - cu' - (K* - u)(l - a — -) 

K* - K*v + v 
= g» ^ - eg' - (1 - 1 - a - + 

1 + kK*(l — w) 

_ u T= fav — JT* + K*v — v + 1 + kK*(l - g) 
= 1 + ^(1-.) a] 

akK*v 

= -(I-V) ; < 0. 

V > \ + kK*{l-v) - 
We then verify the second component in (|1.2p . 

K* -K*v + v , 



cw' + t>[l 



l + fcif*(l-«) J 



-v(l — v) + v[l — 



(l + kK*(l-v))(l-k + ak) l + kK*{l-v)> 

. " . - w(l + fcX* - K*) - " , (1 - u) 

r 1 — k-\-ak v / 1 — fc+a/c V ' i ^ 

[ l + fcl^l-u) J ~ ' 
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since 1 + kK* — K* 



l-k+ak " 

Thus (u,v) forms a smooth upper-solution for (|1.2[) . 



□ 



We next construct the lower solution pair for system (|1.2p . For a small but fixed 
number I with < I < 1 — k + ka, we begin with yet another KPP system: 



v" - cv' + 



(2.8) 



(1 - k + ak)(l + kK*(l - lv)) 



1 + kK* - IK* . 
v(l v) = 0, 



1 + kK* - K* 



1 + kK* - K* 

v( — 00) = 0, V(+00) = ; ; < 1. 

v ' V ' 1 + kK* - IK* 



Corresponding to the notions in Lemma [2] 



/(*) 



(1 - k + ak)(l + kK* (I - lv)) 



1 + kK* - IK* , 
^~ l + kK*-K* ^ >0 



for v e (0, It^lZ^ )- /(0) - /( jjfe^ ) = 0, f (0) = a > 0, and f( jjfe^ ) 

1 — i+ia ^ n 

l-Z+2a(l-fc+afc) v u - 

For each fixed c > 2y / a, define 



(2.9) 



«(0 
«(0 



K*lv{0 

m 



with ?;(£) a solution of (|2.8[l . 



Lemma 4. for eflsc/i c > 2-y/a , 



smooth lower solution of system 11. fy) . 



Proof. On the boundary, one has 



and 



u(-oo) \ / 
v(-oo) / I 



w(+oo) 
w(+oo) 



1 + fciT" -K* 

K I 

1 + kK* - IK* 

1 + kK* - K* 
V 1 + kK* - IK* J 



< 



K* 
1 
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Furthermore, 

" i iv* K*-u + v 

u -cu-(K -u)(l-a- -) 

1 + k{K* — u) 

= K l[v — cv — (- — v)(l — a 



I _/v l + k(K*-u) 

* r , 1 w K* - U + V , 

= if*/[-( T - - a- ~ - 



, l + kK*-lK* 



(1 + kK*(l - lv))(l - k + ak)~ y l + kK*-K 

i+fc(g*-ro 

u[l + kK* - K * - (1 + fei^* - *)v] 
1 + k(K* - u) 



} 



K * 1 v_{{\-v)[K*l-l-k{l-a)K*l]+ a 



1 + fc(if* - u)~ I 1-k + ak 1-fc + afc 

«[1 - k + ak - 1 - T ^ + - kv + k(l - a)lv] 



l + k(K* -u)- r " "" 

> 0, 

l-fc(l-a) 2 



since < Z < 1 — k + ka < * ffp "A and 



v" - cv' + v(l - —-=- — =-r) 

- - v 1 + k(K* - u) ' 



a 1 + kK* — IK* K* — u + v 

-v[l v) + v(l — 



{l + kK*{l-lti)){l-k + ak)- y l + kK*-K*~' - v l + k(K*-u) 

K* — K*lv + v a l + kK*-lK* 

1 + k(K* -u) ~ (1 + kK*(\ -lv)){l-k + ak) ^ 1 + kK* - K* 

l + kK*(l- Iv) -K*+ K*lv~v~ T3 fe(l - \%\ K K '.-J£ v), 



- L l + kK*(l-lv) 

1 + fcX* - K* - -r 

= v[l 



1 + kK* - K* - 1 _ fc % fc - v(l + klK* + K*l - (1 + kK* - LPT*)) , 



l + kK*{l-lv) 



2 *A-*(l-0 



- l + kK*{l-lv) ~ ° 
Therefore, the conclusion of the lemma follows. □ 



We next show, by shifting the upper solution far to the left, that the upper- and 
the lower-solutions as derived in Lemma [3] and Lemma |4] are ordered. 
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Lemma 5. Let c > 2~Ja be fixed and (u,v) T , (u,v) T be respectively the upper and 
lower solutions defined in (|2.7[) and (12. 9p . £/ien there exists a number r > 0, suc/i 
i^a* (i2,v) T (£ + r) > (u,v) T {£) for £ e E. 

Proof. The proof is only for c > since the one for c = 2 v / a is similar. We first 
derive the asymptotic behaviors of the upper- and lower-solutions at infinities. By 
Lemma [2] 

/ u \ ( K*A 1 

(2.10) I j (0 = [ | , ^ + o(, ; ) 

and 





as £-> +oo, where /z = |(c- yc 2 + 4 1 _, + ^ ( ' 1 H k +a fc) ) < and A x , A 2 , S x , B 2 are 
positive constants. 

Since (|2.6|) is translation invariant, v r {£,) = v(£ + f) is also a solution of (|2.6p 
for any f € 1. It then follows that (u r , v r ) T (^) is an upper solution for system 
(|1.2p . For the asymptotic behavior of (u,v) r (£,) at — oo, we can simply replace 
(K*A U Ai) T by (K*A 1 , A{) T e^ '^*-^* in (|2TT0|) . Now we choose f > large 
enough such that 



(K*A U A 1 ) F e^ r > (K*A l7 A x f . 

Then there exists a sufficiently large N\ > such that 

/ u f (0 \ ( u(0 \ 

(2.14) > far£e(-oo,-J\Ti]. 
V / V ««) / 

On the other hand, the boundary conditions of the upper and lower solutions at 
+oo also imply that there exists a number JVjj > such that 

/ u f (0 \ / u(0 \ 

(2.15) > for£e [N 2 , +oo). 

We next show that the inequalities (|2.14j) and (|2 . 15[) also hold on the interval 
[— Ni, N 2 }- There are two possible cases to deal with: 



TRAVELING WAVE SOLUTION 



Case 1. If we already have 

(2.16) > | on [-JVi.JVa], 

we then let f = r and obtain the Lemma. 

Case 2. There exists a point £ e (— JVi, A^) such that 

/ « f (fo) \ / H(&) 
(2-17) < 

V ^(eo) / V «te>) 

with strict inequality holding for at least one of the two components. 

In this case, we will use the the Sliding Domain method. We first shift (u r , v r ) to 
the left by increasing f until we can find a n > f > such that (it ri (£), ?J ri (£)) T > 
(u(£),v(£)) T on the interval [-iVi, 7V 2 - (n - r)]. We then shift (u ri (£), u ri (£)) T 
back to the right by decreasing ri to some r 2 > f such that one of the branches of the 
upper solution touches its counterpart of the lower solution at some point £ 2 in the 
interval (— Ni+r 2 , N% — (r% — f)). On the endpoints of the interval (— iVi + r 2 , iV 2 — 
(fx — r)), we still have (u T2 (£), ■u T ' 2 (£)) T > (u(£), n(£)) T . In summary, we now have 
" r2 (6) = «(&) and iT 2 (£) > u(0, « r2 (0 > «(0 f or £ e (-iVi + r 2 , N 2 - (n - f)). 

Letting W(Q := (u T2 ,v T2 ) T (£,) — (u,v) T (£) and F = (F%, F 2 ) T = (-(if* - u)(l - 
a - T§#l^), «(1 " TS^ry )) T For £ e (-N, + r 2 ,N 2 - (n - f)) we have, 



(2.18) W - cW' 



V7 = 0, 



for some C« G [0,1], « = 1,2,3,4. Since the above system is monotone and the 
cube [(0,0), (K*, 1)] is convex, we can readily deduce by Maximum Principle that 
> for £ e [— iV~i + r 2 ,N 2 — — f)]. This is a contradiction, i.e., such £ 2 
does not exist and we can further decrease r 2 to r. This shows that £o does not 
exist either. We therefore have (u r ,v r ) T '(£) > (it, v) T (£) for £, € K. □ 

To ease the burden of notations, we still use (u, v) T to denote the shifted upper 
solution as given in lemma [5j With such constructed ordered upper and lower 
solution pairs, we now have, 

Theorem 6. For every c > 2y/a, system il.2]) has correspondingly a unigue (up 
to a translation of the origin) traveling wave solution. The traveling wave solution 
is strictly increasing and has the following asymptotic properties: 

1. Corresponding to the wave speed c > 2\fot, 

P.I., (* (f) W A, V« + »(e«) 
V »({) / \A 2 J 

as £ — > — oo; and 

i2 - m LHJ-ur +o(e > 

as £ — > +oo, where A\, A 2 , A\, A 2 are positive constants. 

2. Corresponding to the wave speed c cr j+i ca l = 2y / a, we have 
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(A Uc + Ai 2c £) 
K* \ ( i n 



1 / V A 22 

-oo, w/iere Ai2 C , ^22c < 0, Au c , A 2Xc e R and An,A 22 > 0. 

Proof. Starting from the upper and lower solution pairs obtained in Lemma [3j4] 
and using the monotone iteration scheme given in [4j 126] , we obtain the existence 
of the solution (u(£),v(£)) T to the first two equations in (|1.2[) for every c > 
which satisfies 

24(0 \ / «(0 \ / «(0 

«(0 / V «(0 / V «(0 

On the boundary, it is easy to see that the solution tends to (0,0) T as £ — > — oo, 
and (if*, 1) T as £ — > +oo according to the last inequality. 

To derive the asymptotic decay rate of the traveling wave solutions at ±oo, we 
let c > 2y/a and 

(2.23) [/(£) := («(£), u(£)) T for - oo < £ < oo 

be the corresponding traveling wave solution of (|1.2p generated from the mono- 
tone iteration. We differentiate (|1.2[) with respect to £, and note that ?/'(£) := 
(wi, w 2 ) T (^) satisfies 

(2.24) w'l - cw[ + An (u, v)wx + A 12 (u, v)w 2 = 0, 

(2.25) w 2 - cw' 2 + A 2 i(u, v)wx + A 22 (u, v)w 2 = 0, 
where 

a i \ -, K* — u + v 

Au(u,v) = 1 



{K- 



l + k(K* -u) 

kv — 1 



A 12 (u,v) 
A 21 (u,v) 
A 22 (u,v) 



'(l + k(K* -u)) 2 ' 
K*-u 



l + k(K*-v,y 

kv — 1 
~(l + k{K* - u)) 2 ' 

K* — u + v 
~ 1 + k(K* - u) 

1 



-v. 



l + k{K*-u) 

Lemma [5] implies that the upper- and the lower-solutions as derived in Lemma 
[3] and Lemma @] have the same asymptotic rates at — oo. (|2.19[) and (|2.2ip then 
follow from Lemma [5j 
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We next study the exponential decay rates of the traveling wave solution L/(£) 
at +00. The asymptotic system of (|2.24p and (|2.25[) as £ — > +00 is 



(2.26) 



xp'l — cip[ — atpi = 0, 

^ - afa + (1 - fctyi - ip 2 = 0. 



It is easy to see that system (|2.26[) admits exponential dichotomy [5]. Since the 
traveling wave solution (it(£), v(£)) T converge monotonically to a constant limit as 
£ — > ±00, the derivative of the traveling wave solution satisfies (wi(±oo),w 2 (±oo)) = 
(0, 0) ([26], p658 Lemma 3.2). Hence we are only interested in finding bounded so- 
lutions of (|2.26p at +00. 

Since the first equation in (I2.26|) is decoupled, we can write its general solution 

as 

V>i = A e 2 e + B 1 e 2 5 

for some constants A 1 and B . Since w\ — > as £ — > +00, one immediately has 
A 1 = 0. 

We then study the second equation of (|2.26[) . rewriting the equation as 
(2.27) (fo)" - c(ih)' ~i>2 = -(1 - Qih, 



since c Vc 2 +4a j g no ^. ^ e cnarac teristics of the homogeneous equation, one has the 
following expression for the general solution of (|2.27p . 

(2.28) <^ 2 =B x e a ^ + B 2 e 2 ? + B 3 e 2 £. 

Since w 2 (£) — > as £ — > +00, then B3 = 0. Also noticing that (|2.27|) is non- 
homogeneous, we have B\ ^ 0. By roughness of the exponential dichotomy ([5]) 
and comparing with the upper solution, we obtain the asymptotic decay rate of the 
traveling wave solutions at +00 given in (|2.20p and (|2.23[) . 

We next show the strict monotonicity of the traveling wave solutions, which will 
be a key ingredient in locating the eigenvalues of the linearized operator in the next 
section. By the monotone iteration process [26], the traveling wave solution t/(£) is 
increasing for £ 6 1, it then follows that (w 1 (^),w 2 (0) T = U'(£) > and satisfies 

// , dFi dFi 

(2.29) w 1 — cw 1 + — — (u,v)wi + ——(u,v)w 2 = 0, 

ou ov 

dF 2 dF 2 

(2.30) w 2 - cw' 2 + ——(u,v)wi + -7^—(u,v)w 2 = 0, 

Ou Ov 

and 



(2.31) (wi(0,M0) T > 0, (w llW2 ) T (±oo) = 0. 

The strong Maximum Principle implies that (wi,w 2 ) T > for £ e R. This 
concludes the strict monotonicity of the traveling wave solutions. 

On the uniqueness of the traveling wave solution for every c > 2-y/a, we only prove 
the conclusion for traveling wave solutions with asymptotics (|2.19[) and (|2.20p , since 
other case can be proved similarly. Let ?7i(£) = (ui,vi) T and U 2 (£) — (u 2 ,v 2 ) T be 
two traveling wave solutions of system (jl.ip with the same speed c > 2y/a. There 
exist positive constants Aj, Bi, i = 1,2,3,4 and a large number N > such that 
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for f < -N, 
(2.32) 



/ c + Vc2-4o 

/ A\e 2 



+ o(e ^ 5) 



V A 2 e 



:+\Z<: 2 -4q 



^ 3 e 5 C 



(2.33) C/ 2 (0 
and for £ > N, 

(2.34) C^(0 = 



c+v / c 2-4e 



o e 



Ate 5 5 



/ 



X* - J5ie 



c-V / c 2 +4n , 



o(e 2 e ), 



V 1-B 2 e 

X* - B 3 e 5 1 » 



(2.35) 



u 2 (0 



\ 1 - B 4 e 



:-Vc 2 +4n , 



+ o(e 2 ? ) 



/ 



The traveling wave solutions of system 1|1.2[) are translation invariant, thus for any 
6 > 0, C/i (£) := £/i(£ + 6*) is also a traveling wave solution of (TO)) . By ([232]) and 
(I2.34[) , the solution [/j (£ + 0) has the asymptotics 



(2.36) C7?(0 
for £ < -JV; 

(2.37) £/?(£) 



A 2 e 



- + Vc 2 -4q c + V c 2 — 4o 

2 p e 2 



c+Vc 2 -4a f , 

+ o(e 5 £) 



2 °e 2 



K* -Bxe 



., c 2 +4a a c- % /c 2 +4q f \ 

2 e 2 



V l-B 2 e 



for £ > N. 

Choosing 9 > large enough such that 



(2.38) 
(2.39) 
(2.40) 
(2.41) 

then one has 
(2.42) 



e 2 



> A 3 



+ o(e 2 £ ) 



A 2 e 2 > At, 



Bie 5 6 < B 3 , 



B 2 e * e < £4. 



> ^2(6 
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for £ e (-00, -JV]U[JV + 00). We now consider system (|1.2I) on [— JV, +JV]. First, 
suppose Uf (£) > C/ 2 (£) on [-JV,+JV], then the function VF(£) = (wi(£), w 2 (£)) T := 
U( (£) - C/ 2 (£) > and satisfies for some 6 (0, 1), i = 1,2, 3, 4, 
(2.43) 



W" - cW" 



^(w2 + C3^1,Wl), ("1,^2 + C4W2) 



1^ = 0, £g(-JV,JV), 



W(~N) > 0, W(+JV) > 0. 



Since the above system is monotone, we can readily deduce by the Maximum Prin- 
ciple that W(g) > on [— JV, N]. Consequently, we have C/f(£) > J7 2 (£) on R in 
this case. 

Second, suppose there is some point in (— JV, JV) such that one of the components, 
say the j-th component, satisfies (E/f (£))j < (C^2(0)j at that point, j = 1 or 2. 
We then increase 9, that is shift £/f(£) further left, so that Uf(-N) > U 2 (-N), 
Uf(N) > U 2 {N). By the monotonicity of E/f and Z7 2 , we can find a G (0, 27V) 
such that in the interval (— JV, JV), we have E/f (£ + 0) > E/ 2 (£). Shifting E/f(£ + 0) 
back until one component of E/f (£ + 0) first touches its counterpart of E/ 2 (£) a t some 
point £ G [— JV, JV]. Since E/f (£ + 0) and E/ 2 (£) are strictly increasing, f G (— JV, JV), 
while at £ = ±JV, we still have E/f (£ + 0) > E/ 2 (£). However, by Maximum Principle 
for that component again, we find that component of E/f and U 2 are identically 
equal for all £ G [— JV, JV] for a larger than the original one for which (12.42j) holds. 
This is a contradiction. Therefore, we must have 

uno > m) 

for all £ G M, where is the one chosen by means of (|2.38p ~ (|2.41[) as described 
above. 

Now, decrease until one of the following situations happens. 

1. There exists a 9 > 0, such that E/f (£) = E/ 2 (£). In this case we have finished 
the proof. 

2. There exists a 9 > and £1 G R, such that one of the components of E7 
and E/ 2 are equal there; and for all £ G R, we have E/f (£) > E/ 2 (£). On applying 
the Maximum Principle on R for that component, we find E/f and U 2 must be 
identical on that component. To fix ideas, we suppose that the component is the 
first component. Then E/f — U 2 satisfies (|2.29|) and (I2.31| . Plugging W\ = into 
f)2 . 29(1 again we find that there is at least one £0 such that W 2 (£g) = 0. Then by 
applying maximum principle to (|2.30|) . we have u> 2 (£) = for £ G R. We have then 
returned to case 1. 

Consequently, in either situation, there exists a 9 > 0, such that 

U?(0 = E/ 2 (£). 

for all £ G R. □ 

The next Theorem shows that the lower bound 2y/a for the wave speed c is 
optimal, hence it is the critical minimal wave speed. 

Theorem 7. There is no monotone traveling wave solution of (|1.2p for any < 
c < 2y/a. 

Proof. Suppose there is a constant c with < c < 2y / a and a solution V(£) = 
(wi,f2) T (£) of (|1.2p corresponding to it. Similar to the proof of Theorem |6l the 
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asymptotic behaviors of V(£) at — oo are described by 

\ v 2 (0 J \B S j \B S j 

where (A s , B S ) T and (A S ,B S ) can not be both zero due to the non-homogeneity of 
the limit system of (|2 . 29[) and (|2.30[) at — oo. The condition < c < 2 v / a implies 
that V(£) is oscillating. This concludes that any solution of (II. 2|) with c < 2^/a is 
not strictly monotone. □ 

Remark 8. An interesting implication of Theorem [H] is that the populations u and 
v have the same exponential rate determined by a, the measurement of the penalty 
for altruism that u receives. The bigger the penalty, the faster the population u 
decreases and population v increases. Given the same population growth rate of u 
and v , according to the model, the cooperators always loses the game. 



3. Stability of the Traveling Waves 

In this section, we investigate the stability of the traveling waves with non-critical 
speed. We first show that the traveling wave solution obtained in Theorem [5] with 
non-critical speed is unstable in the continuous function space C(M) x C(M). This 
motivates us to investigate the spectrum of the linearized operator (see (|3 . 3|) ) in 
the exponentially weighted Banach spaces. 

Retaining all the transformations and rescalings in section [2J we consider system 
(|1.2[) in the following equivalent form, with initial conditions as given in 

r u t (x,t) = u xx + f(u), 

\ U(x,0) = U(x). 

where U = (u, v) T and F(U) = (-(K* -u)(l-a- I jggj=^ T ),„(l- 1 ^"+^ )) T , 

U(x) is the initial function. 

Rewriting the above equivalent system in terms of (£, t) variable, with moving 
coordinates £ = x + ct, we have 

( o U / Ut&t) = Ux-cUt + F(U), 

[ j I u(to) = u(0. 

For a fixed wave speed c* > 2y/a and £ = x + c*t , U*(£) = (u*(0> v * (0) T 
obtained in Theorem[6]is then a steady state of (13. Let [/(£, t) = U*(£) + V(£, 
be a solution of (|3.ip . We then obtain the following system for the perturbation 
function V(£, t); 



( . 9 , f v t cv + n(v,u*), 

where 

(3.3) C V = Vtt-c'V e + ^(U*)V 
is a linear operator, and 

dF 

(3.4) K(V, U*) - F(C/* + V) - F(U*) - ^j{U*)V 
is a nonlinear operator. 
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The stability of the traveling wave solution J7*(£) in a specific Banach space 
depends critically on the location of the spectrum o~(£) of C Since C is defined on 
R, it has point spectrum (eigenvalue) 

cr p (£) = {A e cr(£) | A is an isolated eigenvalue of finite multiplicity} 

as well as essential spectrum a e (C) = a(C)\a p (C) ( [2j H3 H3 [21] ) . Let C(M) be 
the space of all continuous functions on K and Co be a subspace of C(M), 



C = {U(0 1 17(0 e C(R) x 



lim ||(7(0l|-0} 
l£|->°° 



where 



U\\ Co = sap \\ U(0 



Theorem 9. The traveling wave solution U*(j£) of system \3.1\) with non-critical 
wave speed c* > , obtained in Theorem^ is unstable with initial functions in 

C . 



Proof. We prove the trivial solution of (|3.2p is unstable. Thus, it suffices to show 
that in the space Co the operator C in (|3.3[) has essential spectrum with positive real 
parts. It turns out that ([7J [53]) the essential spectrum of the operator C is bounded 
by the spectrum of £ at ±oo, which are denoted by C + and C~ respectively. More 
precisely, we have 



OF 

C+V = Vx-c*V i + —(U* + )V 



(3.5) 



V e6 -c*V e + 



1 - k 



v, 



(3.6) 



V^-c*V ( 



Vx-c*V 6 



dF 
dU 



(U1)V 



-(l-a)(l-fc 




ak) 1 — a 



where U± denotes the limit of U* (0 as £ — > ±oo. 
Now consider the equation 

at 

Following [53], we replace V by e^ XT+t ^^I, where I is the identity matrix. We then 
have 



(3.7) 



r) F 

3 (A t +,«)(_ c2/ _ c * cu + ^(U* + ) - XI) = 0. 



The spectrum of the operator C consists of curves 



(3.8) 

Therefore we have 

(3.9) 

or 



BF 

det(~( 2 I - c*£il + ^{UX) -XI) = 



C 2 - c*(i - a - A = 0, 



(3.10) 



C 2 - c*(i - 1 - A = 0. 
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Let A = x + yi, then by (|3.9p we have 
(3-11) x = - T -rz-a, 

or by (l3~TuT) . 

(3.12) x = __^__i. 
Similarly, the spectrum of C + consists of curves: 

(3.13) x = -j^-(l-a){l-k + ak), 
or 

(3.14) x = --f- +a 

in the complex plane. Consequently, we have 

maxi?e a e (C) — max{— a, — 1, —(1 — a)(l — k + afc), a} = a > 0. 
Hence, by [7J [2U] , the steady state U* (£) of (|3.ip is unstable in Co . 

□ 

Theorem|9]says that the Co norm for the initial perturbation is too large to stabilize 
the traveling wave solutions. To have further control on the exponential rates of 
functions in Co, we introduce the following weighted Banach space. Let u\, oi be 
two non-negative numbers and the space Cj 1i(T2 be defined as: 

C„ u a 2 = MO | U(0(e a ^ + e~^) e Co}, 

with the norm 

imK,„ 2 =sup|| U(0(e^ + e-^) || . 

(2) 

We can also similarly define the Banach space C^,^ as: 

ci% 2 ={u\ u(-),u'(-),u"(-)ec auai2 } 

with norm 



|t/|L(2) =supS l=0 



(2) 

It can be easily verified that Cj 1i(J2 and Cri- ljC r 2 are both Banach spaces. 

We now restrict the initial conditions and the operator C to the newly defined 
Banach space C CTliCT2 with ci > 0, cr 2 > and a\ + o\ ^ 0. To relate the operator 
C in Co to an equivalent operator in C ffl (T2 , we introduce the following translation 
operator T : C ai .<j 2 — + Co defined as 

TV := {e ai Z + e - a ^)V. 

The operator T is thus linear, bounded and has a bounded inverse T -1 : Co — > 
C CT1 , CT2 with r- x V = (e CTl « + e- a ^Y x V. 
Consider operator 

(3.15) CV = TCT-W. 
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One can easily see that £ is a linear operator with domain C^ 2 ^(R) x C^(R). By 
relation (|3.15p . considering C in C ai ^ 2 is equivalent to considering C in Co- The 
operator C can be explicitly written as 

dF 

(3.16) CV = V ee - (2 ffl + c*)Vj + (2. 9 2 _ 52 + c * 5l + ^(C/*))V, 

where (2gi + c*) and (2g 2 — 52 + c*Si) in (|3.16p represent the matrices (2gi + c*)I 
and (2gf - g 2 + c*gy)I with 



.91 (0 
.92(0 



-(72^ 



and further properties 



lim 



£ — >oo 



5i(C) = o-i j lim ? ^_oo.9i(0 = -cr 2 ; 



lirru_ 



»52(0 = of, 



lim^-oo^CO = of- 



For the convenience of later use let M(£) = 2c?^ — g 2 + c*gi + §jj(U*). We now 
locate the essential spectrum of the operator C in the space Co- 

Lemma 10. Let c* > 2y / a and if <J\ and er 2 satisfy 



,„ „ s „ -c* + \/c* 2 + 4a c* - Vc* 2 - 4a c* + Vc* 2 - 4a 
(3.17) 0<ai< , <a 2 < , 

then the essential spectrum of the operator C in the space Cq (R) x Co (R) is contained 
in a closed sector in the left half complex plane with vertex on the horizontal axis 
left to the origin. Outside this sector, there are only eigenvalues and the resolvent 
oft. 

Proof. Same as in the proof of Theorem^ we first study the operator C at infinities. 
We have 



(3.18) 
(3.19) 



dF 

C+V = Vx - (2*1 + c*)V ? + (a 2 + c** x + ^-AITW, 



C~V = % - (-2<t 2 + c*)V £ + (a 2 - c** 2 



dU 

dF 
dU 



(U*))V, 



where a\ +c**i + j£j{U*) and of — cV 2 
matrices 



Qfjj{U*) are, respectively, the constant 



(3.20) 

and 

(3.21) 



M+ = 



rj + c *i — a 

1-k 



c*o~i — 1 



M~ 



c**2 — (1 — a)(l — k + ha) 



1 - a 



ct 2 — c* (T 2 + a 

The essential spectrum of the operator C + consists of curves: 
det(-C 2 ^ - i{2*i + c*)IQ + M+) = 0. 

Therefore, 

-C 2 - i(2*i + c*)C + o\ + c**i - a - A = 0, 
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or 

-C 2 - i(2(Ti + c*)C + o\ + c*a ± + e x - 1 - A = 0. 
Equivalently , one has 

or 

y 2 , _2 , „* 



+ c a\ + e\ — 1. 



(2^+c*) 2 

Similarly, the essential spectrum of £~ consists of curves: 

V 2 

x = ~ I o ■ *\2 + a 2 ~ C * CT 2 - (1 - a)(l - * + ka), 



Consequently, we have 

max{Re<7 e (£)} = max{er 2 + c*o\ — a, a\ + (fax — 1, 

(3.22) 

cr| — c*(T2 — (1 — a)(l — + ka), a\ — c 02 + a} 
A simple calculation shows the validity of the Lemma. □ 

Corollary 11. Assume that all the hypotheses in Lemma VHK hold, then the essential 
spectrum of the operator C in the space C ai ,a 2 ^ es * n the left half complex plane. 

Proof. The conclusion follows immediately from Lemma [TOl and relation (|3.15p . 

□ 

We now locate the eigenvalues of the operator C in the space C ai , CT2 . The next 
lemma shows that is not an eigenvalue of the operator C in the space C CTl)CT2 . 
The idea of the proof is as follows. If there is an eigenfunction associated with 
the eigenvalue zero of the operator £ in the space C CTl)0 - 2 , then we can show the 
derivative of the traveling wave solution is larger than or equal to any multiple of 
this function by examing their asymptotics at infinities. But the boundedness of 
the derivative of the traveling wave solution makes this impossible. 

Lemma 12. Assuming all the hypotheses in Lemma \Tb\ then is not an eigenvalue 
of the operator C in C ffl . ff2 . 

Proof. Let U* (£) be a traveling wave solution of (|3.ip found in Theorem [5] By 
taking Gateaux derivative we can verify that 

(3.23) C((U* (0)') = 0. 

Suppose that there exists a nonzero function V £ C rTl . a2 satisfying the equation 

(3.24) CV = 0. 

we then claim that the inequality |rF(£)| < (£/*)'(£) is consequently true for all 
r e K and £ e R. Writing S = {r e R\ \rV(£)\ < (U*)'(£), £ € R}, we verify the 
following properties of S: 

1. S is non-empty since € S. 

2. S is closed. Let € S, i = 1, 2, ... and r, ^ r as i — > +00, then we will have 
\nV(0\ < ( u *Yi which implies that \rV(£)\ < (U*Y(Q, we therefore have r e S. 
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3. S is open. Let r G S, we will show that there exists a S > such that 
(r - 6, r + 6) c 5. We cfatm that |rV"(0| < (^*)'(0 implies |rV"(0| < (^* )'(£)■ 
In fact, let W(£) = (17*)' (0 - then ^(0 > for £ G R and it satisfies the 

following equation: 



(3.25) 



w" — cw[ + Auwi + A\2W2 = 0, 
w 2 ' - cw' 2 + A 2 \wi + A22W2 = 0, 
Oi,w 2 )(-oo) = (W!,W 2 ) (+00) = 



where j4y, i, j = 1,2 are the entries of the Jacobian §jj(U*). The Maximum 
Principle implies that W(£) = {wi(0,w 2 (0) T > for £ G R. This_shows that 
(U*)'(0 - rV(0 > 0, f G R. Similarly we can show that (t/*)'(£) + rV(£) > for 
The claim then follows. 

We next show that the claim further implies |fV"(£)| < (U*)'(£) as long as r is 
sufficiently close to r. According to condition (|3 . 1 T[) and the assumption that V G 
C CTli(73 , for any fixed f G R, there exists a N > sufficiently large such that (e ^ + 
e"^) [(£/*)'(£) - rF(£)] > for all £ G (-00, -N], this implies that (Z7*)'(£) > 
rV(£) also holds there. Furthermore, due to the claim and the boundedness of the 
functions (U*)' and V", we can find a <5 > such that for any f G (— 6 + r, 6 + r), 
([/*)'(() > fV"(£) on the finite interval [— iV, iV]. Now we fix f — f and show 
[(£T)'(0 - ff(0] > for £ G [JV, +00). 

Noting that the entries of the diagonal of the matrix §^(t/*(+oo)) are both 
negative, we can choose a vector P + > such that (by increasing N if necessary) 
i£(t/*(£))P + <0for£G[iV,+oc). 

We have to deal with two cases: 

Case A. If we already have (*/*)'(£) - fF(£) > for £ > iV, then the Maximum 
Principle implies that (£/*)'(£) - rV(£) > on [N, +00). Analogly (*/*)'(£) + 
rV(£) > is also true for £ G R. Consequently, S is open. 

Case B. If there is a point in the interval (—00, — N) such that one of the com- 
ponents of vector (£/*)'(£) — rV(^) takes negative local minimum on this point, 
we consider function W(£) = (f *)'(£) — rV(£) + tP+. The asymptotic rates 
of (If*)' and V at +00 imply that there is a sufficiently large r > such that 
W = (*/*)'(£) - rV(£) + tP + > for £ G [iV,+oo). We further assume that one 
of the components of W(£), say W\ for example, takes minimum at a finite point 
£2 in [N, +00). It is not hard to verify that there is a r = T2 > such that the 
corresponding W(£) staisfying wife) = and W(£) > for £ G (— 00,— N). For 
such t 2 on the one hand, we have 

CW = - c*W ( + ^j{U*)W 

(3.26) 

= r 2 ^(U*)P+ <0, 

on the other hand at £ — £,2, the first component to the left hand side of (|3.2G|) 
is larger than or equals zero. We then have a contradiction, and consequently 
([/*)'(£) — fV(£) > for £ G [N, +00). We are again in the situation descibed by 
case A. By a similar argument, we can show (t/*)'(£) +rV(£) > for £ G [N, +00). 

In summary, for any f G (— 8 + r, 6 + r), we have |fV(£)| < (U*)', i.e., S is open. 

Now the set 5* is a non-empty, open and closed subset of R, then S = R, but this 
is impossible since (U*)' is bounded. Therefore the equation (|3.24p does not have 
non-zero solution in C ai t<72 . □ 
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The next lemma shows that there is no eigenvalue of C in C ai _ a , 2 with positive 
or zero real part. 

Lemma 13. Let Cq be the complexified space of Cq(M.) and A be an eigenvalue of 
the operator £ with U_ £ Cq as its eigenf unction, then Re A < 0. 

Proof. Let A = Ai + \ 2 i and J7(£) = U 1 ^) + iU 2 (£) for (el, where Ai, A 2 S K 
with \\ + Al ^ and W(0 G C* (K). 
Consider the Cauchy problem ([3"1 128j). 

(3.27) F^ry-AiF, V^(e,0) = C/ 1 (e). 

It is easy to verify that V(£, t) = t/ 1 (() cos(A 2 t) — ?7 2 (£) sin(A2i) is bounded and 
solves (|3.27l) for £ S K and t > 0. We may assume that at least one of the 
components of V takes positive value at some point of £ and t, otherwise we can 
consider —V. Since the vector T(U*)'(£) — > +oo as £ — > — oo, similar to the proof 
of Lemma [TJl we can choose a large £o > such that 

(3.28) M (£)P + < 
for £ > £o, and at the same time 

(3.29) V(f > *)<T(IT) , (0 

for £ < -£ and t > 0, and M(£) = 2 3 2 - g 2 + c*5i + §tj(U*) in (|3~2^]) . Moreover, 
we can choose a small I > such that 

(3.30) (e 2 /-e(2si+c*)I + M(£))P + < 0. 

for £ > £o- We can choose an even larger £ if necessary such that at least one 
of the positive values of V occurs inside the interval (— £cb£o)- The positivity of 
T( [/*)'(£) for £ e M implies that there is a r > such that 

(3.31) V&t) < rT(C/*)'(0 fore e [-&,&], * > 0. 

Suitably adjusting r such that the equality in (|3 . 3 1 1) holds on at least one component 
at a point (£i,*i) with £i 6 [— £o>£o] and ti > 0. Now let Ai > 0, we have 
Claim: V(f , t) < rT(U*)' for all £ € R, i > 0. 

The claim is already true for £ < £o by the discussion above. Suppose the 
assertion is not true for £ > £q, then we can find a £ > £cb and a t\ > such that 
> rT(Z7*)'(£). Since Q+(£) = e^P+ -► +oo as £ -► +oo, there is a f > 
such that V(£, t) < rT(J7*)'(£) + fQ+(£) for all £ > £ and t > 0, and for at least 
one j (j = 1,2) and one £2 > £0 and one £2 > 0, we have the equality for the j— th 
component: 

V5(6,*a) = rr(^)'(6) + 

Let F(£,i) = rT(U*)'(£) + fQ + (£) - V(£,t), then Yj has the following properties: 
*i(&,*2) = 0, 19 (&,*) > 0, Y^,*) > for all £ > * > 0. Hence Y jtt ^ 2 ,t 2 ) < 0, 
X?',?(^2,*2) = and Yj^i&ffo) > 0, and t) satisfies 
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Yj,t = ~V ht 

= -(LV-XiV^ 

(3.32) > (-£V + \ 1 V + £rT(U*y + jC.rQ + -\ 1 (rT(U*)' + rQ+)) j 
= (CY~\ 1 Y) 

= Y m - {2g x + c*)Y M + M j: )\ + M jt2 Y 2 - \) r 

owing to (HI30I), CrT(U*Y = and rT{U*)' + fQ + ) > 0. However, we find that 
at (£2,^2) the left hand side of the inequality (|3.32p is smaller than or equal to 
zero, while the right hand side is larger than zero. We then have a contradiction. 
Therefore the claim is proved. 

We next show that the above is also true even for r = 0. In fact, let S = {r > 
0|V(f,i) < rT{U*)'{Cl for £ G K, t > 0}, by the claim S ^ 0, then r = minS > 0. 
We need to show rg — 0. Suppose on the contrary that we have rg > , then 

(3.33) V(f, t) < r r(f7*)'(0 for $ e M, i > 0. 

Since T(t/*)'(£) — > +00 as £ — > -co, the strict inequality holds at £ near -co. We 
first assume that an equality occurs in (|3.33|) at one of the components (say, the 
i-th component) at a point (£, t) with (€l and t > 0. Moreover from 

Y,t > (fy-Aiy)i 

= - (2gi + cf)Y il£ + M i; i(oy + M ii2 (€)y 2 - Ai^ 

with Y(f,t) = r(f7*)'(f) - V(£,i) and (Mi,i, Mi j2 ) the i-th row of the matrix 
function M(£). The positivity Lemma ([23]) implies that Y;(£,£) > for £ S R 
and t > i, further by the periodicity of V, we have that Y (£, t) > for £ S K and 
t > 0. This contradicts with the existence of £ and shows that 

V(& t) < r r(C/*)'(£) for £ e M, i > 0. 
Similar to the justification of the claim in the proof of Lemma [T2l we can show 
there exists a sufficiently small S with r ^> S > 0, such that 

V(f,t) < (r - *)T (ir)'(f) for (61, t > 0. 
It then follows that rg — 5 € S thus contradicts with the definition of rg. Then we 
conclude that rg = but this again is impossible since by assumption a component 
of V is positive at some point £ and t. Hence we have proved the Lemma. □ 

Lemma 14. Under the hypotheses of Lemma \l(K the operator C in Cg is sectorial. 

Proof. See [7]. □ 

Theorem 15. Under the hypotheses of \l(A the operator C generates an analytical 
semigroup in C„ ltCr2 . 

Proof. The conclusion follows from Lemma [14] and [7] . □ 

We now state the stability theorem, 

Theorem 16. Assume the hypotheses of Lemma \l(A The traveling wave solution 
U* of (jl.2p . with wave speed c* > , is asymptotically stable according to norm 
|| • || := || • ||c CT1 „ 2 - That is, if the initial condition [/(£, 0) = J7(£) € C with 
(£/(£) — U*(£)) € C<x li0 - 2 and \\U — U*\\ sufficiently small, then the solution [/(£, t) 
to (|3.1[) exists uniquely for all t > and satisfies 

(3.34) \\U(€,t) -U*{x + ct)\\< Me- bt , 
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where the constants M > 0, b > are independent oft and U. 

Proof. Corollary [11] Lemma [TO] Lemma [12] and [13] show that the spectrum of 
the operator C in the space C ai _ a2 is contained in an angular region in the left 
complex plane. Thus, the trivial solution of system (|3 . 1 [) is linearly asymptotically 
stable. Theorem [T5] further implies the local nonlinear stability of the traveling 
wave solutions ( [2 H3 [20] ) . □ 

Acknowledgement: The authors thank Professor Anthony W. Leung of Uni- 
versity of Cincinnati for the inspiration of the lower solution derived in Lemma 
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